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ABSTRACT 

We present a generalization of the concept of magnification bias for 
gravitationally-lensed quasars, in which the quasars are selected by flux in more 
than one wavelength band. To illustrate the principle, we consider the case of 
two-band selection, in which the fluxes in the two bands are uncorrelated, per- 
fectly correlated, or correlated with scatter. For uncorrelated fluxes, we show 
that the previously-held result — that the bias is the product of the single-band 
biases — is generally false. We demonstrate some important properties of the 
multi-band magnification bias using model luminosity functions inspired by ob- 
served correlations among X-ray, optical, infrared and radio fluxes of quasars. In 
particular, the bias need not be an increasing function of each flux, and the bias 
can be extremely large for non-linear correlations. The latter fact may account 
for the high lensing rates found in some X-ray/optical and infrared/radio selected 
samples. 

Subject headings: Cosmology: gravitational lensing 

1. Introduction 

If a massive galaxy lies along the line of sight to a background quasar, the galaxy may 
act as a gravitational lens, magnifying and forming multiple images of the quasar. Beginning 
with the pioneering work of Turner, Ostriker, & Gott (1984), many authors have computed 
the number of lenses that should appear in well-defined samples of quasars, with particular 
attention given to the dependence of this statistic on the vacuum energy density (Turner 
1990; Kochanek 1996; Helbig et al. 1999; Sarbu, Rusin, & Ma 2001; Li & Ostriker 2002). 
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These calculations must take into account not only the probability that a massive galaxy 
will be aligned closely enough with a background quasar (the lensing cross-section), but 
also the enhancement of the quasar flux due to lensing (the magnification bias). This is 
because quasar samples are usually defined by observed flux in some wavelength band, and 
gravitational lensing boosts the observed flux, thereby sampling a fainter portion of the 
quasar luminosity function. For example, if intrinsically faint quasars are sufficiently more 
numerous than bright quasars, then a quasar with a given observed flux is more likely to be 
lensed than the cross-section alone would imply. 

More recently attention has turned toward the information about galaxy mass profiles 
that can be gleaned from lens statistics. These statistics include lensing rates (see, e.g., 
Keeton & Madau 2001; Wyithe, Turner, & Spergel 2001; Li & Ostriker 2002), the ratio of 
four-image to two-image lenses (see, e.g., Rusin & Tegmark 2001; Finch et al. 2002), the image 
separation distribution (see, e.g., Kochanek & White 2001) and the brightness distribution 
of central images (see, e.g., Rusin & Ma 2001; Keeton 2001, 2002; Evans & Hunter 2002; 
Oguri 2002). All of these applications of lens statistics require a good understanding of 
magnification bias. 

Borgeest, von Linde, & Refsdal (1991) noted that quasar samples selected by both radio 
and optical flux measurements are subject to what they called a "double magnification bias." 
If the radio and optical fluxes from a given quasar are nearly independent, then quasars 
bright in both bands are especially likely to be lensed 3 . By assuming that gravitational 
lensing produces only one possible value of magnification, and using power-law luminosity 
functions for the optical and radio bands, Borgeest, von Linde, & Refsdal (1991) showed 
that the resulting two-band magnification bias is the product of the bias factors computed 
separately for each band. 

It is timely to revisit the issue of multi-band magnification bias with a more general 
approach. With the advent of large-area sky surveys at many wavelengths, it has become 
possible to define samples of thousands of quasars by their observed fluxes in X-ray, optical, 
infrared, and radio bands. Quasars appear in large numbers in, for example, the RASS 
(ROSAT All-Sky Survey: Truemper 1982; Voges et al. 1999) and eventually ChaMP (Chandra 
Multi-wavelength Project: Wilkes et al. 2001; Silverman et al. 2002) at X-ray wavelengths; 
NVSS (NRAO-VLA Sky Survey: Condon et al. 1998) and FIRST (Faint Images of the Radio 
Sky at Twenty centimeters: Becker, White, & Helfand 1995; White et al. 1997) at radio 



3 Note that the important property of the two bands is independence, not a large separation in wavelength 
(as has since been stated in the literature; Bade et al. (1997)), although of course these two properties are 
related. 
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wavelengths; 2MASS (Two Micron All Sky Survey: Kleinmann et al. 1994) at near-infrared 
wavelengths; and SDSS (Sloan Digital Sky Survey: York et al. 2000; Schneider et al. 2002) 
at optical wavelengths. Cross-correlation of these catalogs (see, e.g., McMahon et al. 2001; 
Ivezic et al. 2002) will become an increasingly important source of information about quasars 
in general, and gravitational lens statistics in particular. 

A few lenses have already been discovered using multi-band selection criteria, at lensing 
rates that are larger than the 0.2-1% typical of single-band lens surveys. Bade et al. (1997) 
discovered the gravitational lens RX J0911. 4+0551 by matching RASS sources with optical 
sources from Schmidt plates. Of the ~ 40 radio-quiet X-ray-luminous high-redshift quasars 
known, two are lensed (Wu, Bade, & Beckmann 1999). A search for very red quasars through 
the matching of FIRST and 2MASS has identified two gravitational lenses out of thirteen 
sources (Gregg et al. 2002; Lacy et al. 2002). None of these projects were designed explicitly 
to discover gravitational lenses, although this is a realistic possibility for the future. 

In this paper we investigate the magnification bias for quasar samples defined by mea- 
surements in multiple wavelength bands. After presenting the basic formalism for iV bands 
(§2), we specialize to the case of two bands and consider some illustrative examples. We 
consider the cases in which the two fluxes are uncorrelated (§2.1), perfectly correlated (§2.2), 
and correlated with non-zero scatter (§2.3). We then use a realistic model of the optical lu- 
minosity function for quasars to demonstrate a few interesting properties of the multi-band 
magnification bias (§3); in particular, the bias does not necessarily increase with flux in 
each band, and there is a profound difference between the case of a linear correlation and 
a non-linear correlation with flux in another band. Finally, in §4 we summarize our results, 
and discuss possible applications of this formalism to real quasar samples. 



2. Magnification Bias and the Multiple Imaging Rate 

We begin by reviewing the case of single-band magnification bias (Turner 1980; Turner, 
Ostriker, & Gott 1984). In a sample of quasars at redshift z with (apparent 4 ) luminosity Li, 
the fraction of multiple-image lensed quasars is 

F(L u z)^- (1) 

where r mu it is the cross-section for multiple imaging, and Bi(L\) is the magnification bias. 
For T mu it <C 1 and -Bir mu i t <C 1, this reduces to the usual expression F(Li, z) = BiT nm \ t . The 



4 By "apparent," we mean that Li is the luminosity inferred from the observed flux and the luminosity 
distance, without taking into account the possible magnification due to lensing. 
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magnification bias is evaluated as 



Bl(Ll '" ) = ' (2) 



where $i(Li, z) is the quasar luminosity function, fi is the sum of the unsigned magnifications 

dp 



of the multiple images, and 4^ is the probability distribution for /i, taken for a singular 



isothermal sphere throughout the paper (dP = 8fi~ 3 dfi for fi > 2). This expression can be 
understood as a likelihood ratio. The denominator is the likelihood that the quasar is drawn 
from the sample of unlensed quasars with luminosity L\ (within dL\). The numerator is the 
likelihood that the quasar is drawn from the fainter sample of quasars with luminosity L ± / fi 
(within dLi/fj,), summed over all possible values of fi. 

Understood this way, the generalization to N bands is straightforward. We require 
knowledge of the multivariate luminosity function, $n(£i, L 2 , L 3 , Ln, z). For a point 
source, the magnification is the same for all bands, because gravitational lensing is achro- 
matic 5 . The multi-band magnification bias is therefore 

„ f^dn^ r ^ 1 ... N (L 1 /n,L 2 /n,L 3 /n,...,L N /n,z) 

±>1..JV(,M, L 2 , L 3 , L N , Z) - — r . [6) 

*1...AT^1, ^2, 1^2,1 J^N, Z) 

The dependence of B on the apparent luminosities of the quasars depends on the corre- 
lations, if any, between the intrinsic luminosities of the quasars in those bands. To illustrate 
the interesting properties that can result, in the following sections we concentrate on the 
simplest non-trivial case, the two-band magnification bias. All of the results are easily gen- 
eralized to N bands. 



2.1. Two-Band Magnification Bias: No Correlation 

Borgeest, von Linde, & Refsdal (1991) considered quasars observed at both optical 
and radio wavelengths, and assumed a power-law luminosity function for each band. They 



5 Gravitational lensing magnification is sensitive to source size. Therefore if the emission regions for 
the two bands differ greatly in their spatial extent, then there is the possibility of wavelength-dependent 
magnification. For example, the small optical emission region of a quasar may be microlensed by stars in the 
lensing galaxy, whereas the more extended emission regions at infrared or radio wavelengths is not generally 
microlensed (see Wyithe & Turner (2002) for a recent discussion of microlensing and magnification bias). We 
ignore the possibility of microlensing in this paper, but note that since the mean magnification of microlensed 
sources equals the magnification of un-microlensed sources, the results presented will be qualitatively correct, 
even for correlations involving bands that are subject to microlensing. 
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showed that if the fluxes in these bands are statistically independent, and if there is only 
one possible value of the lensing magnification, then the two-band magnification bias is 
equal to the product of the biases that would be computed separately for the optical and 
radio bands. This result is not true in general. As we show below, even if the two bands are 
independent, the result does not hold because real gravitational lenses produce a distribution 
of magnifications. 

First, we reproduce the result of Borgeest, von Linde, & Refsdal (1991) using our for- 
malism. For N = 2, Eq. (3) is 

B 12 {L U L 2 , z) - — — — - — . (4J 

5*12(^1, ^2, Z) 



If the bands are independent, then $ 12 (Li, L 2 , z) = $i(Li, z)Q 2 (L 2 , z). The lens model used 

dP 



by Borgeest, von Linde, & Refsdal (1991) can be described by ^ = 5(fi — /i ), in which case 



rj It t x_ 1 $l 2 (L 1 //i ,L 2 //X ,z) _ 1 $i(Li/// ,z)$ 2 (L 2 ///0,^) _ , . , s 

(5) 

This results fails for the more realistic case in which there is a range of possible mag- 
nifications, because ^ appears once in the numerator of the multi-band magnification bias, 
but appears separately in each numerator in the product of the single-band biases. For 
example, following Borgeest, von Linde, & Refsdal (1991), suppose $i(Li,z) = ^i^L" 1 and 
$ 2 (L 2 , z) = $2^X2 2 . If we adopt the magnification distribution appropriate for an isothermal 
sphere = for all /i > 2), then 

for cti + a 2 > —4. Analogous calculations of the single-band bias factors give 

B±(Li, z)B 2 (L 2 , z) = \ 6 2 -(4 +ai+a2); (?) 

which can be either larger or smaller than Eq. (6). 



2.2. Two-Band Magnification Bias: Perfect Correlation 



If the two bands are perfectly correlated, with L 2 = f(L\), one might expect that 
no new information is provided by the observation in the second band, and therefore that 
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the two-band bias is equal to the single-band bias for either band. This is not quite true. 
Gravitational magnification multiplies both fluxes by the same factor. If the unmagnified 
fluxes are linearly correlated, then the magnified fluxes also obey the correlation. However, 
if the correlation is non-linear, then the magnified fluxes do not obey the correlation, and the 
source must be gravitationally lensed. In the the appendix, we derive this result formally, 
by calculating the magnification bias for general correlations (see the next section) in the 
limit of zero scatter. 



2.3. Two-Band Magnification Bias: Imperfect Correlation 

More generally, L\ and L 2 are correlated with some intrinsic scatter. On physical 
grounds we expect the magnitude of the scatter to scale with the luminosity [i.e. AL2/L2 ~ 
g(ALi/Li)], which makes it convenient to use logarithmic variables I = logL. Suppose that 
the correlation between L x and L 2 is a power-law, L\ = Li, or l x = 7/2- Because of the 
correlation, it is convenient to express the luminosity function in terms of the new variables 
U\ = Hi + I2 and U2 = —7/1 + h, which describe the location parallel and perpendicular 
to the correlation, respectively (see Fig. 1). In these variables, the luminosity function 
(expressed in density per square logarithmic interval) is ^12(^1, u 2 , z) — (7 + ^)$i2(h, h, z). 
The luminosity function can also be written ty 12 (ui,U2, z) = ^i(ui, z)p(u2\ui, z), where 
^i(ui, z) is the luminosity function in the new variable U\ and p{u2\u\, z) is the conditional 
probability of U2 given u\. Because we expect the scatter to be symmetric in reflection about 
the correlation 6 , we assume p(v,2\ui,z) to be Gaussian with variance a, hence 

*i2(ui,«2,*) = #i(ui,z)-^=^exp (-^)- ( 8 ) 

Defining M = log //, the magnification bias is 

dP <S> 12 (h-M,l2-M,z) 



B 12 (h,l2,z) = / dM 



oc 



dM $12(h,l2,z) 

[ 1 \ ■ ^ - L (9) 



dP (7 + ^)*i2 k - (1 + $)M, u 2 + (7 - 1)M, z] 

dM 



dM (~f+±)y 12 (u 1 ,U2,z) 



6 As an example of why we expect the scatter to be symmetric in reflection about the correlation, consider 
a sample of quasars with flux measured in two optical wave bands, say r and i. We would expect to find a 
variation of bias with i-band at fixed r-band that is qualitatively similar to the variation with r-band at fixed 
i-band. This symmetry in the magnification bias requires symmetry of the scatter in reflection about the 
correlation. We therefore choose to model the scatter as a symmetric function in logarithms of luminosity; 
that is, defined normal to the correlation (i.e. along the u-i axis). 
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Inserting Eq. (8), 

n , s r° „w dP *iK -M',z] ( 1 \f 7(7-l)»w\ 2 2 

" (10) 

where we have defined M' = (1 + ^)M. This expression illustrates many important points. 
First, if the correlation is linear (7 = 1) then B\ 2 is independent of u 2 , and the contours of 
constant bias run normal to the correlation (i.e. along lines of constant l 2 + h). Furthermore, 
if z) is a monotonically decreasing function of ui, then we find that B 12 {li, h, z ) is an 

increasing function of both li and l 2 (since u 2 increases monotonically with both li and l 2 ). 
This example will be further explored in Case 1 of §3. Eq. (10) also demonstrates the behavior 
arising from non-linear correlations (7 7^ 1). Here the exponential plays an important role; 
it introduces an asymmetry in the bias across the correlation. If, for example, 7 > 1, then 
large biases can result from negative values of u 2 (i.e. below the correlation), because the 
exponent becomes positive. On the other hand, the exponent is negative for all u 2 > 0, and 
hence the bias above the correlation is small. This example will be further explored in Case 
3 of §3. 

In the following section, we evaluate Eq. 10 numerically, with more realistic assumptions, 
in order to illustrate these and other interesting and potentially observable properties of the 
multi-band magnification bias. 



3. Magnification Bias for Illustrative Bi-Variate Luminosity Functions 

We consider measurements made in two bands, and a power-law correlation between 
the two bands: L 2 = Lj. As in §2.3, the scatter (normal to the correlation) is assumed to be 
Gaussian in logarithms with half- width a. We consider 4 examples of Eq. (8). The first two 
examples involve linear correlations (7 = 1), and the second two examples involve non-linear 
correlations. 

For the first of the two correlated bands, we use a luminosity function $i(Li, z) that is 
appropriate for optical wavelengths. A good representation of the observed optical quasar 
luminosity function at redshifts z < 3 is provided by the following double power-law form 
(Boyle, Shanks, & Peterson 1988; Pei 1995): 

rh (T \ ®*/L*(z) 

° ( ' Z) ~ [L/L.{z)]* + [L/L m (z)]*>- (U) 

At the faint end, the logarithmic slope of this function is —f3\ = —1.58, while at the bright 
end the slope is — (3 h = —3.43 (Boyle et al. 2000). Moreover, all dependence on redshift (for 
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z < 3) is in the break luminosity L*(z). We therefore show the luminosity function in units 
of L* throughout the remainder of this paper. Setting $i(Li) = $ (Lx), we find from Eq. 8 
that & (Li) is related to ^i(u 1: u 2 ) through 



This equation defines the functions \E , i(m 1 ) used in this section. 

1.7 = 1.0, a = 0.15: A linear correlation with constant scatter. This is the situation that 
might be expected between two different optical bands. The contours of $12 (gray), 
and of B12 (black), are shown in the top left panel of Fig. 2. As was derived in the 
previous section (and might be expected intuitively), the bias increases monotonically 
with both luminosities, and is constant along lines normal to the correlation. 

The top right panel shows the corresponding single-band luminosity function (gray 
line), and magnification bias (dotted black line). These functions are the same for 
both bands because of the linear correlation. In addition we plot _B 12 (L 1 ,L 2 ) for two 
paths through (Li, L2)-space: one at fixed L\ (thick dashed line, plotted as a function 
of L 2 ) and the other below but parallel to the correlation (thin dashed line, plotted as 
a function of L x ). 

2. 7 = 1.0, <y = 0.15 — 0.02 log -ui: Same as the previous example, except in this case we 
allow the logarithmic scatter to depend on luminosity. The results are shown in the 
lower panels of Fig. 2, in the same format as the previous example. 

The contours of magnification bias wrap around the contours of $12, increasing rapidly 
as one moves normal to the correlation (along the u 2 axis). This can be understood as 
follows. Magnification draws quasars from regions of lower intrinsic luminosity, where 
the probability density of quasars (as given by $12) is larger. This is especially true 
when the observed luminosities fall at some distance from the correlation, because the 
scatter in the correlation is larger at lower luminosities. 

The dependence of B\ 2 on L 2 , for fixed Li, is again shown by the thick dark dashed line. 
Interestingly, the dependence is not monotonic. For small values of L 2 the magnification 
bias is large. The bias decreases as L 2 rises through the expected intrinsic value, 
and then increases again. The bias along the path denoted by the thin dashed line 
demonstrates that the bias can become very large for sources below the correlation. 

3. 7 = 1.5, o = 0.2: A non-linear correlation with a constant scatter. This situation 
approximates the correlation that has been observed between the X-ray (L 2 ) and optical 




(12) 
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(Li) bands for quasars (Brinkmann et al. 2000). Results for this case are shown in the 
top two panels of Fig. 3, in the same format as the previous examples. 

In this case, $2(^2) (thin gray line) is a flatter function than <3>i(Li) (thick gray line), 
and therefore the single-band bias B 2 (thin dotted line) is smaller than Bi (thick dotted 
line). Although B12 is an increasing function of L\, it is actually a decreasing function 
of L2 (for fixed L\). This runs counter to the naive expectation that the brighter the 
quasar is (regardless of band), the more likely it is to be lensed. The reason is that 
when L 2 is smaller than expected from the correlation, reducing both luminosities by 
the same factor \i (along a line of unit slope, in the top left panel) brings one to a 
region of much higher probability. 

4. 7 = —1.0, a = 0.3: An anti-correlation with constant scatter. This is a somewhat 
artificial example, but we might imagine there are two ways for a quasar with a fixed 
energy source to radiate its energy, and one of these ways can be blocked by a variable 
amount. For example, the optical and far- infrared luminosities might be expected to 
exhibit some degree of anti-correlation due to dust obscuration. 

The results are plotted in the lower panels of Fig. 3. The contours of magnification 
bias are parallel to the anti-correlation. For small luminosities the bias is smaller than 
unity. Quasars in this region are less likely to be lensed than the cross-section alone 
would imply, because lensed quasars would be drawn from a population with very small 
density. Conversely, for large luminosities, the bias becomes arbitrarily large. 

4. Discussion 

The multi-band magnification bias is an a posteriori statistic. It is used to estimate the 
probability that the apparent luminosities of a given quasar, as measured in several bands, 
are due to gravitational magnification, rather than being intrinsic to the quasar. When a 
sample of quasars is selected through the matching of sources in two different catalogs, both 
fluxes must be used to perform this calculation. One must also have some knowledge of the 
intrinsic correlation (if any) of the fluxes, and the distribution of magnifications produced 
by lensing. 

One might expect that the multi-band bias is maximized when the bands are uncor- 
rected (an example of which is shown for the radio-optical correlation of SDSS early data 
release quasars in Fig. 4), since in that case there is no redundant information in the flux 
measurements. Upon further reflection, or using the mathematics developed in this paper, 
one realizes that this is not true — the relevant information is how discrepant the observed 



-10- 



fluxes are from the correlation, and whether the discrepancy can be made smaller if the 
observed fluxes are all reduced by a constant factor. 

Many of the illustrative examples presented in this paper approximate certain correla- 
tions that have been observed for real quasars. In particular, the multi-band magnification 
bias may result in very high lens fractions for certain quasar samples. First, we consider the 
case of a quasar sample selected by optical colors. The top left panel of Fig. 4 is a logarithmic 
plot of SDSS 2-band vs. r-band fluxes for the SDSS early data release quasars (Schneider et 
al. 2002). The data show a linear correlation with scatter, the magnification bias for which 
is illustrated in the top panels of Fig. 2. The magnification bias must be computed using 
both optical measurements, unless the sample is 100% complete in one filter (i.e., unless 
after selecting quasars in i, the r-band magnitude was measured in every single case). As 
an example consider the sample of SDSS z > 5.8 quasars (Fan et al. 2001). Since the z- 
band selection is at ~ HOOA in the rest-frame, quasars with fixed absolute B magnitude 
(~ 4400A) are more likely to be selected if they are bluer than average. Thus a sample of 
quasars selected in this manner will be bluer than average and lie blueward of the correlation 
on a plot of the intrinsic correlation between My and M-q. The magnification bias for these 
sources may be significantly smaller than that computed using only extrapolations of the 
i?-band luminosity function (Wyithe & Loeb 2002; Commerford, Haiman, & Schaye 2002). 

Next, we consider examples of non-linear luminosity correlations. Brinkmann et al. 
(2000) measured the correlation between ROSAT X-ray and FIRST radio fluxes for matched 
quasar samples. They find that while radio-quiet quasars show a linear relationship between 
X-ray and radio luminosity, radio-loud quasars have an X-ray flux L x that varies with the 
radio luminosity L T as L x oc L® A8±om with an intrinsic scatter of ~ 0.2 dex. Furthermore, 
Brinkmann et al. (2000) showed that the X-ray luminosity correlates non-linearly with optical 
luminosity L , following L x oc £,i- 42 ± - 09 . The second of these correlations (in flux) is plotted 
in Fig. 4 for quasars in the SDSS Early Data Release (Schneider et al. 2002), but can be 
seen more clearly in Fig. 14 of Brinkmann et al. (2000). The multi-band magnification 
bias corresponding to the second correlation 7 is illustrated in the top panels of Fig. 3. The 
magnification bias can be extremely large for sources that are luminous at both optical 
and X-ray wavelengths. This may be the explanation for the apparently high probability 
of lensing in bright X-ray selected quasar catalogs (Bade et al. 1997). The location of the 
gravitational lens RX J091 1.4+0551, which was selected from cross-correlation of optical and 
X-ray catalogs, is shown on this plot by the large dot 8 . The fluxes place the quasar below 



7 Note that while we have presented results for a non- linear correlation with index 7 > 1, the result for 
7 < 1 is simply obtained through reversal of the axes. 

8 We used the integrated R from Bade et al. (1997) and color transformations from Fukugita, Shimasaku 
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the correlation, in the region where we expect the magnification bias to be large (see the 
upper panels of Fig. 3). The lens HE 1104-1805 is also X-ray loud (Wisotzki et al. 1993; 
Reimers et al. 1995). While this lensed quasar was not discovered through cross-correlation 
between catalogs, it is interesting to note its location on this plot, shown by the open square 
in Fig. 4. The quasar is found to be very bright in both bands, and is again in the region of 
high magnification bias. It is suggestive that the two X-ray loud gravitational lenses both 
appear to lie in the region of high magnification bias, as expected. 

Finally, the multi-band magnification bias may also provide an explanation for the large 
gravitational lens fraction (2 out of 13) found through the matching of FIRST and 2MASS 
sources (Gregg et al. 2002; Lacy et al. 2002). Fig. 4 shows the correlation for near infrared 
luminosities verses radio luminosities compiled from table 1 of Barkhouse & Hall (2001). The 
radio/near-IR correlation appears to be steeper than linear. If true we might expect very 
large biases for luminous near-IR sources. The top panel of Fig. 3 demonstrates magnification 
bias for a non-linear correlation, and shows that the bias of around 100 necessary to achieve 
a lens fraction of 2/13 is possible. 

5. Summary 

This paper has discussed the multi-band magnification bias for gravitational lensing 
with arbitrary luminosity functions in several bands. Previous discussion of the multi-band 
magnification bias (Borgeest, von Linde, & Refsdal 1991) focused on the case where the fluxes 
in the two bands are independent. If a single value for the lens magnification is considered, 
they showed that this assumption leads to a multiple magnification bias that is equal to the 
product of the single-band biases. However, we have shown that this equality breaks down 
in the more realistic case when there is a distribution of possible magnifications. 

We also discussed the multi-band magnification bias when the fluxes in the various bands 
are correlated. In the case of a perfect (i.e. zero scatter) linear correlation, the information 
from the second band does not change the magnification bias. However, if the correlation 
is non-linear, then sources with fluxes that obey the correlation cannot be lensed. On the 
other hand, sources with fluxes that do not obey the correlation must be lensed. 

Of course, real correlations have intrinsic scatter. We have calculated the multi-band 
magnification bias for bi-variate luminosity functions with finite scatter about both linear 
and non-linear correlations. For a linear correlation (as expected for a quasar sample selected 



& Ichikawa (1995). 
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by optical colors) we find that the magnification bias is an increasing function of either flux. 
Calculations of lens statistics from incomplete color-selected quasar samples should therefore 
account for the multi-band magnification bias. 

Non-linear correlations (and anti-correlations) with finite scatter were also explored. If 
the fluxes in two bands are correlated through a relation that is steeper than linear, then 
sources that lie below the correlation can be subject to a very large bias. The observed 
correlation between X-ray and optical flux (and possibly between infrared and radio flux) 
for quasars is steeper than linear. Suggestively, the two known X-ray loud gravitationally 
lensed quasars lie below the X-ray/optical correlation in the region of large magnification 
bias. Thus the multiple magnification bias may provide an explanation for the large lensing 
rates found in X-ray/optical and infrared/radio selected samples. 
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A. Two-Band Magnification Bias for Perfect Correlations 



In this appendix we derive the results mentioned in § 2.2 for the two-band magnification 
bias in the case of perfect correlations by taking the limit of small a in Eq. (10). It is 
convenient to rewrite the exponential as follows 



lim B 12 (h,l 2 ,z) 



cr^O 







",M-^^- M '^l im exp 
dM' V^u^z) v^o 



la 2 7 + 1 



7(7 " l) M> + 2u 2 



7 + 1 



(Al) 



First, consider the case of a linear correlation (7 = 1). In this case we find the expo- 
nential function is unity for all a > 0, and Eq. (Al) reduces to 
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Note that in this case we must have l\ = l 2 . We can relate ^i(ui, z) to $i(Zi, z) in the limit 
of small a: 
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Thus, for a perfect linear correlation, we find that 
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which is the single band magnification bias. Therefore, if two luminosities obey a perfect 
linear correlation, the magnification bias is simply equal to the single-band magnification 
bias computed from either luminosity. 

Next consider the case of a non- linear correlation (7 7^ 1). Specifically, we assume 7 > 1 
and < ML- < M' < 00. Then 



lim oxp I l -M'^- ^ 

^0 1 \ 2a 2 7 + I 



7(7 -y +2^ 
7+1 




1 

00 



if u 2 > - 
if u 2 = - 
if u 2 < 



7(7 - 1) M' 
7 + 1 2~ 

7(7 - 1) M' 
7 + 1 2~ 
7(7 - 1) M' 



7 + 



7 +l J "m 

that lim^o B 12 (l ± , l 2 , z) = 00 if w 2 < — 7 ^ 7 ~^ iVf+ n . Sources that lie on the correlation have 



Note that only sources with u 2 = or u 2 < — 7l7 +1 ±j M+ n are allowed. As a result we find 
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U2 = 0, and therefore linv^o Bi 2 (h, h, z) = 0. Thus if the correlation is non-linear, sources 
that lie on the correlation cannot be lensed, while sources that lie off the correlation must 
be lensed. 
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Fig. 1. — Schematic of variables defining the luminosity function for non-zero scatter about 
a power-law correlation. The grey lines are contours of a bi-variate luminosity function and 
are shown for context. 
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Fig. 2. — Bi-variate magnification biases for linear correlations with scatter. Top: Linear 
correlation, with a scatter that is insensitive to luminosity (7 = 1.0, a = 0.15). Bottom: 
Linear correlation with a scatter that decreases with luminosity (7 = 1.0, a = 0.15 — 
0.02 log Ui). The left hand panels show contours of the bi-variate luminosity function (grey 
lines). The solid lines are contours of magnification bias. The right hand panels show the 
corresponding single band luminosity functions (grey lines). Also shown are the single band 
magnification biases (dotted lines), and the magnification bias along the paths denoted by 
the dashed lines in the left hand figure. The bias for the path denoted by the thin dashed 
line is plotted as a function of L 1: while the bias along the thick dashed line is plotted as a 
function of L 2 . Because of the linear correlation, the single band luminosity functions and 
magnification biases are identical for the two bands. 
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Fig. 3. — Bi-variate magnification biases for non-linear correlations with scatter. Top: Non- 
linear correlation, with a logarithmic scatter that is insensitive to luminosity (7 = 1.5, 
a = 0.2). Bottom: Non-linear anti-correlation, with a logarithmic scatter that is insensitive 
to luminosity (7 = —1.0, a = 0.3). The left hand panels show contours of the bi-variate 
luminosity function (grey lines). The solid lines are contours of magnification bias. The 
right hand panels show the corresponding single band luminosity functions (grey lines) and 
the single band magnification biases (dotted lines). Thick and thin lines denote quantities 
in L\ and L 2 respectively. Also shown are the magnification biases along the paths denoted 
by the dashed lines in the left hand figure. The bias for the path denoted by the thin dashed 
line is plotted as a function of Li, while the bias along the thick dashed line is plotted as a 
function of L 2 . 
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Fig. 4. — Correlations in different bands. Top Left: SDSS i-band vs. SDSS r-band flux. 
(Schneider et al. 2002) Top Right: FIRST radio flux vs. SDSS r-band flux (Schneider 
et al. 2002). Lower Left: ROSAT X-ray counts vs. SDSS r-band flux (Schneider et 
al. 2002) for quasars with redshifts larger than 0.5. The large dot in this panel repre- 
sents RX J0911. 4+0551, while the open square shows the location of HE 1104-1805. Bottom 
Right: Radio vs. Near IR luminosity (Barkhouse & Hall 2001). In the left hand panels the 
observed correlation lines are drawn to guide the eye. 



